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Abstract 

We show that a basic quantum white noise process formally reproduces 
quantum stochastic calculus when the appropriate normal / chronological 
orderings are prescribed. By normal ordering techniques for integral equa- 
tions and a generalization of the Araki- Woods representation, we derive 
the master and random Heisenberg equations for an arbitrary Gaussian 
state: this includes thermal and squeezed states. 



1 Quantum White Noises 

It is possible to develop a formal theory of quantum white noises which nev- 
ertheless provides a powerful insight into quantum stochastic processes. We 
present the "bare bones" of the theory: the hope is that the structure will be 
more apparent without the mathematical gore and viscera. 

Essentially, we need a Hilbert space Ho to describe a system of interest. 
We postulate a family of pseudo-operators {a^ ■ t > called creation noises; 
a formally adjoint family {a^ : i > 0} called annihilation noises and a vector 
called the vacuum vector. The noises describe the environment and are assumed 
to act trivially on the observables of Hq . 

The first structural relations we need is 



a- VI/ = 



(QWNl) 



which implies that the annihilator noises annihilate the vacuum. The second 
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relations are given by the commutation relations 




(QWN2) 



Otherwise the creation noises all commute amongst themselves, as do the anni- 
hilation noises. Here k ~ ^7 + is a complex number with 7 > 0. We have 
introduced the functional kernels 6± having the action 



/oo 
/(s)<5± (t-s) 
-00 



(1) 



for any Riemann integrable function /. 

If the right-hand side of (QWN2) was just ^6 {t — s) , then we would have 
sufficient instructions to deal with integrals of the noises wrt. Schwartz func- 
tions; however introducing the (5±-functions we can formally consider inte- 
grals wrt. piecewise Schwartz functions as we now have a rule for what to 
do at discontinuities. In particular, we can consider integrals over simpliccs 
{t > ti > ■ ■ ■ > tn > 0} . The objective is to use the commutation rule (QWN2) 
to convert integral expressions involving the postulated noises to equiva- 
lent normal ordered expressions, so that we only encounter the integrals of the 
5±-functions against Riemann integrable functions. 

We remark that the commutation relations (QWN2) arose from consid- 
erations of Markovian limits of field operators af (A) in the Heisenberg pic- 
ture satisfying relations of the type [a^ (A) , a+ (A)] = Kt-s (A) 6 (t — s) + 
Kt-s (A)* 9{s — t) where 9 is the Heaviside function and the right-hand side 
is a Feynman propagator. 



1.1 Quantum Stochastic Calculus 
1.1.1 Fundamental Stochastic Processes 

For real square- integrable functions f = f {t), we define the following four fields 

A'' (/) := / [at] ' [a-] ' f {s) ds, i,j e {0,1} (2) 
Jo 

where on the right-hand side the superscript denotes a power, that is [a]° = 1, 
[o]^ = a. 

With l[o,t] denoting the characteristic function of the time interval [0,t], we 
define the four fundamental processes 

Ai^ := A'^ (l[o,t]) = /* [aty [a-]' ds, i,j e {0, 1} ; (3) 
Jo 



2 



1.1.2 Testing Vectors 



We denote by the Hilbert space spanned by the symmterized vectors 

A® . • . <E)fn ■■= (/i) ...A^° (/„) The Hilbert spaces are then nat- 
urally orthogonal for different n and their direct sum is the (Bose) Fock space 
^ = ®^=o-^^"^- The exponential vector with test function / e (M+) is 
defined to be 

'if)--=En^^^f- (4) 

Note that (e (/) |e {g)) = cxp7 and that vE- = e (0). 

For a subset T of (K^) , the subset of Fock space generated by the 
elements of T is €(T) = {e (/) : / e T}. In general, we shall denote the 
Fock space over a one-particle Hilbert space f) by F ([)) = ^) i thus 

jr = r(L2 (IR+)). 

1.1.3 Quantum Stochastic Processes 

Let W be a fixed Hilbert space, we wish to consider operators on the tensor 

product Ti ® T. A family of operators (^t)(>o defined on a common domain 
D ® *t (T) can be understood as a mapping from : DxDxTxTx M+ — » C 
:((/), ^, /, 5, t) ^ ((^ (g) e (/) |Xt V £ {g))- 

1.1.4 Quantum Stochastic Integrals 

Let (^t"') be four adapted processes. The quantum stochastic integral hav- 
ing these processes as integrands is 

Xt = f ds [aty Xi^ [a^y (5) 
Jo 

where we use the Einstein convention that repeated indices are summed (over 
0,1). We shall use the differential notation dXt = [a^y Xl'' [a^y dt or even 

i^^[atyxna^y. 

The key feature is that the noises appear in normal ordered form in the 
differentials. Suppose that the i^t^^ defined on domain D (g) 

then it follows that 



'P®e{f)\^i>®e {g)) = [/* m {<i> ® £ (/) ^ ® 6 {g)) [g {t)y 

for all (j),^ € D and f,g € TZ. (The equivalence can be understood here as being 
almost everywhere.) 
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1.1.5 Quantum Ito's Formula 

Let Xt and Yt be quantum stochastic integrals, then the product XtYt may be 
brought to normal order using (QWN2). In differential terms we may write this 
as 

d{XtYt) = {dXt)Yt + Xt{dYt) 

= (dXt) Yt + Xt (dYt) + (dXt) (dYt) (6) 

where the Ito differentials are defined as {dX^ Y, := [a+X Xl^ {Yt) [a^]^' dt; 

Xt (dYt) := [at]\Xt) Y^^ [a,"]' dt; (dX^ (dY^ := \4X X^Y^^ [a,"]' dt. 

The quantum Ito "table" corresponds to following relation for the funda- 
mental processes: 

(dJ^f){dA't>)=^{dA^i). (7) 

1.1.6 Quantum Stochastic Differential Equations 

The differential equation = [a^j'X/"' \a^X '^'ifli initial condition X^ = Xq 
(a bounded operator in Ho) corresponds to the differential equation system 

^®e{f)\^i;^e (5) J) = [f* ® e (/) \Xi^ ^ s {g)) [g {t)Y 

and in general one can show the existence and uniqueness of solution as a 
family of operators on Ho J^- The solution can be written as Xt = xq + 
ds [aff Xl^ [djy ■ In Hudson-Parthasarathy notation this is written as Xt = 
Xo+ /o Xl^i^dAl'' , where the tensor product sign indicates the continuous tensor 
product decomposition = T (i^(0, t]) ® T {L'^{t, 00)) . 

1.2 Quantum Stochastic Evolutions 

A quantum stochastic evolution is a family ( Jt)j>g mapping from the bounded 
observables on Ho to the bounded observables on Ho d) ^- We are particularly 
interested in those taking the form 

Jt (x) = ulxut 

where Ut is a unitary, adapted process satisfying some linear qsde (a stochastic 
Schrodinger equation). In the rest of this section we establish a Wick's theorem 
for working with such processes. 

1.2.1 Normal-Ordered QSDE 

Let Vt be the solution to the qsde 

^=Li,[atXVt[at]\ Vb = l; (8) 
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where the are bounded operators in T-Iq. The associated integral equation 
is Vf = 1 + /q dti Lij [af^y Vt^ [it]]'' which can be iterated to give the formal 



series 



Vt = 1 + ^ / dti...dtn Li^ji ■ ■ ■ Li^j„ 

n=l "'t>tl>---tn>0 

X [<]"•••[<]" K"]^" (9) 
= Nexpj^ dsLij[aty[a;yY (10) 

where N is the normal ordering operation for the noise symbols af. 
Necessary and sufficient conditions for unitary of Vt are that 

i,,+4+7iUy = 0. (11) 

(Necessity is immediate from the isometry condition d(u}Ut^ = U}d{Ut) + 
d (c//) Ut + d (c/t) d{Ut) = Ul (iy + l]^ + 7iLiij) Ut ® dAi^ = 0, but suf- 
fices to establish co-isometry d (UtU^^ = 0.) Equivalently, we can require that 

Ln = -{W-1), Lw = L, 
^ 1 

where W is unitary, H is self-adjoint and L is arbitrary. 

1.2.2 Time-Ordered QSDE 

Let Ut be the solution to the qsde 

^ = ^E,,[aty[a-yU,, C/o = 1; (12) 

where the Eij are bounded operators in TIq. Here wc naturally interpret Tt = 
Eij [afy [a^y as a stochastic Hamiltonian. (For Tj to be Hermitian we would 
need En and Eqq to be self-adjoint while s|q = Eqi.) 
Iterating the associated integral equation leads to 

Ut = l + ^(-z)"/ dh...dtnE,,j,---E,^j^ 

Jt>ti>...tn>0 

x[<]"M" (13) 
= f exp I -z ^ ds Ei^ [at] " [a;] ' | , (14) 

where T is the time ordering operation for the noise symbols af. 
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1.2.3 Conversion BVom Time- Ordered to Normal- Ordered Forms 

Using the commutation relations (QWN2) we can put the time-ordered expres- 
sions in (12) to normal order. The most efficient way of doing this is as follows: 



[at ,Ut] = at ds Eij [a+]' [a^ Y Us 

I Jo 

= -i [ dsEij[at'at] [a-f 
Jo 

= -inEij [a^yUt 
or a^Ut — Utttf = —iKEiio^Ut — inEioUt- This implies the rewriting rule 

aiUt = (1 + inEii)-^ {Utat - inEwUt} . (15) 

ThxisiE,, [a+y [a^y Uc = iE.o [a+y Ut+iEa [a+Y {1 + iKEn)~^ {Uta^ - iKEwUt}. 
From this we deduce the following result: 



Theorem 1 Time- ordered and normal- ordered forms are related as 



Texp 



ds Eij [a+] * [a^ ] ^ = N exp <^ / ds Lij [a+] ' [a^ ] 



(16) 



where 



Lu = -iEii {1 -\- inEu) \ iio = -i (1 + iK-Bn) ^ Ew, 
Loi = -iEqi (1 + iKEii)~^ , Loo = -iEoo + kEqi (1 + iKEu)~^ Ew 



(17) 



2 Gaussian States 

The above construction requires the existence of a vacuum state ^P; by applica- 
tion of creation fields we should be able to reconstruct a Hilbert-Fock space for 
which ^ is cyclic. But what about non- vacuum states? We describe now the 
trick we shall use in order to consider more general states for the simple case of 
one bosonic degree of freedom. 

Let a, satisfy the commutation relations [a, a^] =1. A state ( ) is said to 
be Gaussian or quasi-free if we have 

(exp {iz*a -|- iza^}) = exp |-n (n -I- 1) zz* -\- m* z^ + mz*^ + iz*a -\- zza*| ; 

(18) 

in particular, (a) = a, (aa^) = n + 1, {aa) = m. Prom the observation that 
((a + Aat)^( a + Ao^) ^ > it follows that the restriction |m|^<n(n+l) must 
apply. 
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Now suppose that ai,a\ and 02,03 are commuting pairs of Boss variables 
and let 

o = xai + yal + za2 + a (19) 

where x,y,z,a are complex numbers. The commutation relations are main- 
tained if |a;p — C|j/p + |zp = 1. Taking the vacuum state for both variables 
ai,al {i = 1,2) then we can reconstruct the state ( ) if + \z\'^ = n + 1 and 
yz = m. That is 



n+l-^-^,y = V^, z=^. 20 
n vn 

2.1 Generalized Araki- Woods Construction 

Let2 [) be a fixed one-particle Hilbert space and let j be a anti-linear conjugation 
on f), that is {j(j)\jtp)f^ = for all (l),tp £ f). Denote by A {(j)) the annihilator 

on r (f)) with test function </) (that is, A {(p) e (ip) = ('/>|V')f| ^ (V') )• A state ( ) 
on r (()) is said to be (mean-zero) Gaussian / quasi-free if there is a positive 
operator iV > 0; an operator M with |Mp < AT (AT + 1) and [N, M] = 0; and a 
fixed anti-linear conjugation j such that 

(exp {iA {<!>) + iA^ (</.)}) = exp (0| (2iV + 1) <P)^ - \ {Mcf>\j<P)^ - \ {j<f>\M<t>\ 

(21) 

for all ^ e ^. In particular, we have the expectations 

(A(0)At(v.)) = {mi^)t, 

(A(0)A(V)) = (M</.|.?-^)^ 
(At(<^)^t(^)) = (j^lMV-)^ 

The state is said to be gauge-invariant when M = 0. The case where 
A'' = (1 — e~^^) and M = yields the familiar thermal state at inverse 
temperature (5 for the non-interacting Bose gas with second quantization of H 
as Hamiltonian. The vacuum state is, of course A'^ = 0, M = 0. 

The standard procedure for treating thermal states of the interacting Bose 
gas is to represent the canonical commutations (CCR) algebra on the tensor 
product r (f)) ig) r (f)) and realize state as a double Fock vacuum state. This 
approach generalizes to the problem at hand. For clarity we label the Fock 
spaces with subscripts 1 and 2 and consider the morphism from the CCR algebra 
over r (()) to that over T (^)j ® T {i))^ induced by 

A {(j)) ^ Ai {Xcj)) ®l2 + li®A\ {jYcj)) + li (ZV) (22) 

where 



/ IMP ^ M 

Here we write Ai {tp) for annihilators on F ([)) ■ {i = 1, 2). 



(23) 
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3 Gaussian Noise 

Now let af (t) and (t) be independent (commuting) copies of quantum white 
noises with respective vacua and 'I'2. Let x. y.x be as in (18) and let a be 
arbitrary complex. We consider the quantum white noise (s) defined by 

:= xa^ (t) + ya^ (t) + za^ (t) + a. (24) 
We consider the stochastic dynamics generated by the formal Hamiltonian 

Tt = Cat + C^at + F (25) 
where C and self-adjoint F are operators on Ho- We are required to "normal 
order" the unitary process Ut = Texp|— i /prfsTsj for the given state and 
this means normal order the (t) and the (t). By using the same technique 
as in (13) we find that, for instance, [a^ (t) , Ut] = —i /J ds [a^ (t) , Tg] Ug = 
—iKxCUt- We can deduce the rewriting rules 

{t) Ut = Ut {a^ (t) — iKxC) ; 
a- (t) Ut = Ut (02 {t) - iKZ*C - iKyC^) . 

The qsde ^ = —iTtUt can be then be rewritten as 

^ = -i: TtUt : -iCyUt {-iKZ*C - inyC^) 

-iCxUt i-iKxC) - iChUt ( —inz* C - iKyC^) 

where :TtUt : is the reordering of TfUt placing all creators af {t) and (t) to 
the left and all annihilators aj~ (t) and 03 (*) to the right. Rearranging gives 

^ = ~iC{xa+{t)Ut + yUta^it)+z*a+{t)Ut + a*Ut) 

-iC^ {xUta^ (t) + ya+ (t) Ut + zUta^ (t) + aUt) 

-[iF + K {{n + 1) C^C + nCC'^ + m*CC + mC+Ct)] U^. (26) 

To obtain an equivalent Hudson-Parthasarathy qsde, we introduce quantum 
Brownian motions At = /q (xa{ (s) + ya^ (s) + za^ (s)) ds with the under stand- 
ing that, for Rt adapted, Rt ® dAt = {xRta^ [t] + ya2 (t) Rt + zRta2 {t)) and 
Rt ® dA\ = [xat (i) Rt + yRta2 (t) + z*a^ (t) Rt) . Then 

dUt = -iCUt O dAl - iC'^Ut dAt - GUt dt (27) 

where G = i (F + a*C + aCt) + k ((n + 1) CtC nCCt +m*CC + mCtCt) . 
Note that the quantum Ito table will be 

dAtdA\ = 7 (n + 1) dt; dA\dAt = jndt; 

dAldAl = 'ym*dt; dAtdAt = ^mdt. (28) 
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It is readily shown, cither by normal ordering or by means of the quantum 
stochastic calculus, that the stochastic Heisenbcrg equation is then 

dJt {X) = -iJt ( [X, Ct] ) ® dAl - iJt {[X, C]) ® dAt + Jt {L {X)) ® dt (29) 

where 

L (X) = 7 {(n + 1) C^XC + nCXC^ + mCXC + m*C^XC^} - XG - G^X. 

(30) 

Finally the master equation is obtained by duality: = L' (g) where 

tr{gL {X)}=tr{L'{g)X}. 
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